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Chapter  1 


Introduction  and  Summary  of 
Results 


Near-field  antenna  measurements  have  been  formulated  and  used  extensively  during  the  past 
three  decades  to  obtain  near  and  far  fields  of  antennas  from  near-field  measurements  [1].  In 
particular,  probe-corrected  planar  near-field  meMurements  were  first  rigorously  formulated 
in  1963  by  Kerns  [2],  [3]  using  the  plane-wave  scattering- matrix  theory  for  antennas  and 
antenna-antenna  interactions.  Kerns  derived  the  probe-corrected  formulas  that  give  the  field 
of  an  antenna  (called  the  test  antenna)  everywhere  in  a  half  space  in  terms  of  near-field  data. 
The  near-field  data  is  the  output  of  a  probe,  with  a  known  receiving  characteristic,  obtained 
by  scanning  on  a  plane  in  front  of  the  test  antenna.  The  formulas  take  into  account  the 
receiving  characteristic  of  the  probe  and  assume  that  multiple  interactions  between  the  probe 
and  the  test  antenna  are  negligible.  All  this  work  was  performed  in  the  frequency-domain  so 
that  Kerns’  formulas  determine  the  fields  at  one  frequency  at  a  time.  For  antennas  excited 
by  short  pulses  with  wide  bandwidths  it  would  therefore  be  convenient  to  have  formulas  that 
determine  the  time-domain  field  directly  without  having  to  perform  the  calculations  for  one 
frequency  at  a  time. 

The  purpose  of  this  report  is  to  derive  probe-corrected  planar  near-field  formulas  in  the 
time  domain,  so  that  a  single  set  of  near-field  data  in  the  time  domain  yields  the  field 
of  the  test  antenna  directly  in  the  time  domain.  Formulas  are  derived  for  both  acoustic 
and  electromagnetic  fields  and  the  space  outside  the  region  occupied  by  the  test  antenna  is 
assumed  to  be  isotropic,  homogeneous,  and  lossless. 

The  time-domain  probe- corrected  formulas  are  derived  in  two  different  ways.  First,  they 
are  derived  by  taking  the  inverse  Fourier  transform  of  the  corresponding  frequency-domain 
formulas.  Secondly,  they  are  derived  directly  in  the  time  domain  using  a  time- domain  expan¬ 
sion  of  the  field  and  a  time-domain  receiving  characteristic  of  the  probe.  A  previous  report 
[4]  derived  non-probe-corrected  planar  near-field  formulas  in  the  time  domain  and  numerous 
results  from  that  report  will  be  used  here.  For  example,  the  frequency-domain  plane-wave 
spectrum  formulas  of  [4,  ch.2]  and  the  time-domain  Radon  transform  formulas  of  [4,  ch.3] 
will  be  used  extensively. 

Rccvivad  for  Ptiiltcation  15  April  199A 
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This  report  is  organized  as  follows.  Chapter  2  derives  Kerns’  frequency-domain,  probe- 
corrected  formulas  that  are  relevant  for  the  translation  into  the  time  domain.  Acoustic  fields 
are  treated  in  Section  2.1  while  electromagnetic  fields  are  treated  in  Section  2.2.  (In  the 
electromagnetic  case,  it  proves  convenient  to  define  a  vector  output  of  the  probe.)  The 
derivations  are  performed  by  first  expanding  the  field  of  the  test  antenna  in  terms  of  plane 
waves  [4,  ch.2].  Then  the  definition  of  the  receiving  characteristic  of  the  probe  and  the 
plane-wave  expansion  are  used  to  express  the  output  of  the  probe  as  a  Fourier  integral  of 
the  receiving  characteristic  of  the  probe  and  the  plane-wave  spectrum  of  the  test  antenna. 
Finally,  the  inverse  Fourier  transform  is  used  to  express  the  plane-wave  spectrum  of  the 
test  antenna  in  terms  of  the  Fourier  transform  of  the  output  of  the  probe  and  the  receiving 
characteristic  of  the  probe.  The  special  case  where  the  probe  is  reciprocal  is  also  discussed. 

In  Chapter  3,  we  derive  the  time-domain  probe-corrected  formulas  for  both  acoustic  and 
electromagnetic  fields.  In  Section  3.1,  these  formulas  are  derived  for  the  acoustic  field,  first 
by  using  the  inverse  Fourier  transform  and  the  corresponding  frequency-domain  formulas 
of  Chapter  2.  Then  these  time-domain  probe-corrected  formulas  are  derived  directly  in 
the  time  domain  using  a  time-domain  plane-wave  expansion  of  the  field  radiated  by  the  test 
antenna  and  a  time-domain  receiving  characteristic  of  the  probe.  These  general  time-domain 
probe-corrected  formulas  involve  a  double  spatial  integral  (over  the  scan  plane)  and  a  single 
time-convolution  integral  (over  all  times).  The  integrand  in  the  resulting  triple  integral  is  a 
product  of  the  inverse  receiving  characteristic  of  the  probe  and  the  output  of  the  probe. 

The  special  case  is  considered  when  the  output  of  the  probe,  due  to  an  incoming  time- 
domain  plane  wave,  is  proportional  to  the  time  derivative  of  the  plane  wave.  It  is  shown  that 
for  this  type  of  probe  the  probe-corrected  formulas  simplify  significantly  because  it  is  possible 
to  calculate  the  time-convolution  integral  analytically.  Specifically,  the  far-field  pattern  is 
given  by  a  double  spatial  integral  (which  is  simply  the  Radon  transform  of  the  output  of  the 
probe)  divided  by  the  angular  dependence  of  the  receiving  characteristic.  Finally,  reciprocal 
acoustic  probes  are  considered  in  Section  3.1.  Time-domain  reciprocity  relations  are  derived 
for  these  probes  and  their  time-domain  receiving  characteristics  are  related  to  their  far  fields. 

In  Section  3.2,  the  time-domain  probe-corrected  formulas  are  derived  for  the  electromag¬ 
netic  field.  Also  for  the  electromagnetic  field,  it  is  found  that  the  probe-corrected  formulas, 
in  general,  involve  the  calculation  of  a  double  spatial  integral  and  a  time-convolution  integral. 
For  probes  whose  outputs,  due  to  an  incoming  time-domain  plane  wave,  are  proportional 
to  the  time  derivative  of  the  plane  wave,  it  is  found  that  the  time-convolution  integral  can 
be  calculated  anal3d;ically.  The  far  field  of  the  test  antenna  for  this  type  of  probe  (called 
a  D-dot  probe)  is  given  in  terms  of  a  double  spatial  integral  only.  Finally,  in  Section  3.2, 
time-domain  reciprocity  relations  for  reciprocal  electromagnetic  probes  are  derived. 

Section  3.3  outlines  two  computation  schemes  to  numerically  calculate  the  time-domain 
far-field  patterns  from  sampled  near-field  data  obtained  by  using  the  special  type  of  probe 
mentioned  above.  The  output  of  such  a  probe,  due  to  an  incoining  time-domain  plane  wave, 
is  proportional  to  the  time  derivative  of  the  field  of  that  plane  wave.  The  two  computation 
schemes  are  probe-corrected  analogs  to  the  non-probe-corrected  schemes  presented  in  [4, 
ch.4]. 
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The  first  scheme,  called  the  frequency-domain  computation  scheme,  is  based  on  the 
frequency-domain  formulation  in  Chapter  2.  This  scheme  consists  of  the  following  three 
steps:  (1)  use  the  Fourier  transform  to  calculate  the  frequency-domain  output  of  the  probe 
from  the  time-domain  output,  (2)  use  the  frequency-domain  probe- corrected  formulas  of 
Chapter  2  to  calculate  the  frequency-domain  far  field  from  the  frequency-domain  output  of 
the  probe,  and  (3)  use  the  inverse  Fourier  transform  to  calculate  the  time-domain  far  field 
from  the  frequency-domain  far  field.  This  scheme  makes  use  of  frequency-domain  sampling 
theorems  and  the  fast  Fourier  transform  (FFT). 

The  second  scheme,  called  the  time-domain  computation  scheme,  is  based  on  the  time- 
domain  formulation  in  Chapter  3.  This  scheme  simply  uses  the  formula  that  directly  gives 
the  time-domain  far  field  in  terms  of  the  time-domain  output  of  the  probe  and  its  time- 
domain  inverse  receiving  characteristic.  The  time-domain  sampling  theorem  from  [4,  sec.4.2] 
shows  how  small  the  sample  spacing  between  points  on  the  scan  plane  must  be  to  calculate 
the  far  field  accurately. 

The  conclusions  of  the  comparison  of  two  non-probe-corrected  computation  schemes 
given  in  [4,  sec.4.4]  are  also  valid  for  the  probe- corrected  schemes  of  this  report.  Specifi¬ 
cally,  the  direct  time-domain  computation  scheme  is  much  simpler  to  program  and  use  than 
the  frequency-domain  computation  scheme.  However,  because  the  frequency-domain  com¬ 
putation  scheme  uses  the  FFT  it  is  much  faster  for  large  antennas  than  the  time-domain 
computation  scheme  when  the  full  far  field  is  calculated  for  all  times. 

When  only  part  of  the  far  field  is  calculated,  the  difference  in  computer  time  for  the 
two  computation  schemes  becomes  smaller  and  the  time-domain  computation  scheme  there¬ 
fore  becomes  more  advantageous  because  of  its  simplicity.  Furthermore,  the  duration  of  the 
far-field  pattern  is  extended  erroneously  because  of  the  finite  size  of  the  scan  plane,  and 
this  longer  time-duration  has  to  be  taken  into  account  in  the  frequency-domain  calculation 
scheme.  The  frequency  spacing  in  the  frequency-domain  computation  scheme  has  to  be  cho¬ 
sen  small  enough  so  that  significant  time-domain  aliasing  is  avoided  in  the  calculation  of  the 
time-domain  far-field  pattern.  The  problem  of  choosing  the  frequency  spacing  small  enough 
does  not  occur  for  the  time-domain  computation  scheme  because  no  frequency  spacing  is 
used. 

Furthermore,  the  time-domain  computation  scheme  has  the  capability  to  calculate  the 
far-field  pattern  at  early  times  from  near-field  measurements  taken  at  early  times  only.  This 
capability  is  not  possessed  by  the  frequency-domain  computation  scheme  because  the  near 
field  is  required  for  its  entire  duration  to  calculate  its  Fourier  transform.  For  many  antennas 
fed  by  short  pulses,  the  time  dependence  of  both  the  near  field  and  far  field  consists  of  an 
early-time  part,  which  contains  most  of  the  power,  and  a  late-time  part  which  is  oscillatory 
and  contains  little  power.  The  duration  of  the  early-time  part  may  be  much  smaller  than  the 
duration  of  the  entire  field.  If  only  the  early- time  part  of  the  far  field  is  of  interest,  one  can  use 
the  time-domain  computation  scheme  to  determine  this  part  from  near-field  measurements 
taken  for  early  times  only.  Thereby  one  can  significantly  reduce  the  number  of  near-field  time 
samples  needed  for  the  far-field  calculation.  If  instead  the  frequency-domain  computation 
scheme  is  used,  the  number  of  near-field  time  samples  cannot  be  reduced  because  time 


samples  taken  over  the  entire  duration  of  the  near  field  are  needed  to  calculate  the  Fourier 
transform  of  the  near  field. 

No  matter  which  scheme  is  chosen,  planar  time-domain  near-field  antenna  measurements 
can  eliminate  the  error  in  the  far-field  pattern  due  to  the  finite  scan  plane  because  this  error 
is  separated  in  time  from  the  correct  far-field  pattern.  This  makes  it  possible  to  use  planar 
scanning  in  the  time  domain  to  compute  the  far-field  of  broadbeam  antennas  in  both  the 
time  and  frequency  domains. 

Finally,  a  numerical  example  illustrates  the  use  of  the  time-domain  computation  scheme 
and  one  of  the  acoustic  probe-corrected  formulas.  In  this  example  the  far-field  pattern  of 
an  acoustic  point  source  with  Gaussian  time  dependence  is  calculated  from  near-field  data 
obtained  by  a  nonideal  probe.  The  output  of  this  probe,  due  to  an  incoming  time-domain 
plane  wave,  is  equal  to  the  time  derivative  of  the  plane  wave  times  cos0,  where  9  determines 
the  propagation  direction  of  the  plane  wave.  It  is  found  that  the  time-domain  computation 
scheme  computes  the  far-field  pattern  accurately  and  that  the  result  obtained  by  neglecting 
probe  correction  is  very  inaccurate  for  off-axis  angles  of  observation. 
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Chapter  2 

Frequency-Domain  Formulas 


Probe- corrected  planar  near-field  frequency-domain  formulas  will  be  derived  in  this  chapter 
for  both  acoustic  and  electromagnetic  fields  using  a  simple,  straightforward  approach.  Those 
frequency-domain  formulas  especially  useful  for  the  translation  of  near-field  techniques  to  the 
time  domain  will  be  emphasized.  Throughout  the  report  e"’"*  time  dependence  is  suppressed 
in  all  the  time-harmonic  equations. 


2.1  Acoustic  Fields 

The  excess  pressure  of  the  acoustic  field  in  a  homogeneous  isotropic  stationary  fluid  is  denoted 
by  which  satisfies  the  homogeneous  Helmholtz  equation  =  0  in  every  source- 

free  region.  Here  k  =  u/c  is  the  propagation  constant  and  c  is  the  acoustic  speed.  In 
the  following  the  excess  pressure  $  is  simply  called  the  acoustic  field.  To  use  the  same 
terminology  as  for  electromagnetic  fields,  the  acoustic  source,  which  generates  the  field  $, 
will  be  called  the  test  antenna. 

Consider  Figure  2.1,  showing  the  geometry  for  measuring  the  acoustic  field  $  radiated 
by  an  acoustic  test  antenna  located  in  the  half  space  z  <  zq  and  fed  by  a  single  propagating 
mode  of  amplitude  oq  at  the  reference  plane  5o.  Alternatively,  one  can  think  of  the  test 
antenna  as  a  source  region  that  includes  the  antenna,  the  feed  system,  and  the  generator, 
without  specifying  a  waveguide  mode  with  amplitude  oq.  This  approach  was  taken  in  [4]. 

When  the  test  antenna  is  the  only  object  in  all  space  it  radiates  the  field  given  by  the 
plane-wave  expansion  [3],  [4,  eqs.(2.8),(2.9)],  [5] 

$(f)  =  ^  z>zo  (2.1) 

Zir  7—00  7—00 

where  To  is  the  transmitting  spectrum  of  the  test  antenna  given  by 

A+OO  ^-^-OO 

ky)  =  - -  /  /  +  ^ 

ZTTCZq  7—00  7—00 


(2,2) 


with 

I  ki  +  ki<p 

\  +  k;-  n  ki  +  ki  >  e 

and  f  =  zi  +  yy  +  zi.  The  propagation  vector  for  each  plane  wave  can  be  denoted  by 
k  =  kxX  +  k^y  +  7z.  (The  product  oqTq  equals  the  spectrum  T  defined  in  [4].) 

The  far  field  of  the  test  antenna,  when  it  is  the  only  object  in  all  space,  is  given  by  the 
well  known  expression  (4,  eq.(2.48)] 

gifcr 

^{f)  ~ —iaokcos6——To{kcos<l>s\n0,ksm<f>s'\nO),  z  >  zq  (2.4) 

where  the  spherical  coordinates  {r,6,<(>)  are  determined  such  that  x  =  rcos<^sin^,  y  = 
rsin<^sin0,  and  z  =  rcos0.  It  is  convenient  to  write  the  far  field  in  terms  of  the  far-field 
pattern  fo{0,<f>)  which  is  defined  such  that  the  far  field  is  $  ~  aoJ^o{0,<f>)e''‘’' /r.  The  far- 
field  expression  (2.4)  gives  the  following  relation  between  the  spectrum  To  and  the  far-field 
pattern 

^oi0,4>)  = —ik  cos  0To{k  cos  4) sm0,ks\n<i>s\n0),  0  <  ^  <  7r/2.  (2.5) 

The  goal  of  this  section  is  to  calculate  the  spectrum  To  for  the  fixed  test  antenna  by 
measuring  its  field  on  the  plane  z  =  Zq  using  a  probe  as  shown  in  Figure  2.1.  The  location 
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of  some  reference  point  of  the  probe  is  denoted  by  fg  =  Xoi  +  j/oy  +  and  in  the  following 
the  probe  is  simply  said  to  be  located  at  fo  when  its  reference  point  is  located  at  tq.  For 
planar  scanning  the  probe  translates  without  rotation  in  the  zq  plane. 

Start  by  assuming  that  the  distance  between  the  probe  and  the  test  antenna  is  always 
sufficiently  large  that  multiple  interactions  can  be  neglected.  That  is,  the  reflection  in  the 
test  antenna,  of  the  field  radiated  by  the  induced  sources  on  the  probe,  does  not  afl^ect 
the  output  of  the  probe.  The  effect  of  the  multiple  interactions  can  formally  be  taken 
into  account  by  using  the  plane-wave  scattering-matrix  theory  [3].  However,  this  theory 
does  not  provide  quantitative  information  about  these  multiple  interactions  and  to  obtain 
useful  probe-corrected  formulas  they  have  to  be  neglected.  Therefore,  we  will  not  use  the 
comprehensive  plane-wave  scattering-matrix  formulation  here,  but  instead  use  a  more  direct 
approach  that  exhibits  the  physics  of  the  transmit-receive  process.  Also,  this  direct  approach 
can  be  extended  to  the  time-domain  analysis  in  the  next  chapter. 

The  output  of  the  probe  at  the  reference  plane  Sp  is  denoted  by  6p(fo)  and  the  receiving 
characteristic  of  the  probe  Rp{kj;,ky)  is  defined  to  be  the  output  of  the  probe,  when  it  is 
located  at  fo  =  0  with  a  fixed  orientation  in  the  incident  plane  wave 
This  means  that  the  output  of  the  probe  (when  it  is  located  at  fo)  due  to  the  incident 
plane  wave  is  given  by  Note  that  the  plane- 

wave  field  propagates  in  the  direction  given  by  the  spherical  angles  {0,(f>) 

determined  from  kx  =  A:  cos  ^  sin  d  and  ky  =  ksm<f>sm6  where  k  is  the  propagation  constant. 
When  kl  k^  <  P ,  IS  real  and  the  plane  wave  is  propagating;  when  kl  +  k^  >  7  is 

imaginary  and  the  plane  wave  is  decaying. 

Equation  (2.1)  shows  that  the  plane- wave  component  of  $  propagating  in  the  direction 
(A:*,/:,,)  is  given  by  aQ{2ir)~^To{kx,ky)e*^'‘^^'*''‘»^'*^^^dkxdky  such  that  the  contribution  to  the 
output  of  the  probe  from  this  plane-wave  component  is 

dbp{fo)  =  aoRpikx,  ky)To{kx,  (2.6) 


provided  multiple  interactions  between  the  probe  and  test  antenna  are  negligible.  Also,  we 
are  assuming  throughout  that  the  probe  output  line  is  perfectly  matched,  so  that  there  is  no 
reflection  from  its  termination.  The  total  output  of  the  probe  is  thus  seen  from  (2.1)  to  be 
the  integral  over  the  entire  kx  —  ky  plane  of  the  expression  in  (2.6),  that  is 

/^OO  /'^OO 

/  «p(fcx,  ky)To{kx,  (2.7) 

•00  J^OO 

which  is  the  final  expression  for  the  output  of  the  probe  located  at  fo.  This  plane-wave 
“transmission  formula”  along  with  the  Fourier  transform  formulas  show  directly  that  the 
spectrum  To  for  the  acoustic  field  is  given  by  [3],  [5] 


=  (27)€^:.M  C  C  (2.8) 


That  is,  the  spectrum  To  for  the  acoustic  field  $  is  simply  the  Fourier  transform  of  the 
output  of  the  probe  divided  by  the  receiving  characteristic  Rp  of  the  probe.  This  is  the 
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general  probe-corrected  formula  which  gives  the  spectrum  To  for  the  test  antenna  in  terms 
of  the  receiving  characteristic  and  output  of  the  probe. 


2.1.1  Formulas  for  Reciprocal  Probes 

If  the  probe  is  a  reciprocal  electroacoustic  transducer  [6],  fed  by  an  electromagnetic  waveguide 
in  which  only  a  single  mode  is  propagating,  the  receiving  characteristic  is  given  in  terms  of 
the  spectrum  Tp  for  the  probe  by  the  equation  [6,  eq.(lOa)] 

Rpi^x,  ky)  =  Tpi-kj:,  -ky)  (2.9) 

Voi*>Po 

where  70  is  the  characteristic  admittance  for  the  propagating  mode  in  the  waveguide  feed  and 
po  is  the  static  pressure  of  the  fluid.  Inserting  this  expression  for  Rp  into  the  probe-corrected 
formula  (2.8)  one  finds  that  the  spectrum  for  the  test  antenna  is  given  by 


=  i2.)Z7AZ-k.)  C  C  (2.1D) 


which  is  valid  when  the  probe  is  a  reciprocal  electroacoustic  transducer. 

Finally  let  us  express  the  far-field  pattern  (2.5)  of  the  test  antenna  in  terms  of  the  far- 
field  pattern  of  the  probe.  The  far-field  pattern  of  the  probe  (with  the  reference 

point  of  the  probe  at  fo  =  0  and  the  spherical  angles  6  and  <i>  measured  with  respect  to  the 
{x,y,z)  coordinate  system  of  Figure  2.1)  is  defined  by  $p  ~  apTp[6.,4>)e^^'' (r  where  $p  is 
the  far  field  of  the  probe  and  Op  is  the  amplitude  of  the  probe-input  signal  when  the  probe 
is  used  as  a  transmit  antenna.  Since  the  probe  radiates  into  the  half  space  z  <  Zq,  the 
relation  between  its  far-field  pattern  Tp  and  its  spectrum  Tp  is  given  by  [4,  sec.2.2.1],  [5, 
eq.(18)]  Tp{9^<i>)  =  tfccos^Tp(Acos0sin<>,A:sin^sin^)  for  7r/2  <  B  <  ■k.  This  shows  that 
Tp(— Acos^sin^,  —  A:sin^sin^)  =  [t^cos(7r  — ^)]~*.Fp(7r  — 7r-t-<^),  for  0  <  B  <  7r/2  and  thus 
from  (2.10)  and  (2.5)  the  expression  for  the  far-field  pattern  of  the  test  antenna  is 


ro(B,<^)  = 


TjoPoUlk  cos  B 


/+00  r+00  _ 

/  6p(ro)e~‘ ’■■’■•’dxodyo-  (2.11) 

•00  •/— 00 


(27r)2aoJ’p(7r  *  0,7r  +  0)  y-oo 


Here  f  =  x  cos  <l>sinB  +  y  sin  ^ sin  0  -I-  z cos B  and  Tp{B,  <f>)  is  the  far-field  pattern  of  the  probe 
when  it  is  located  at  fo  =  0  of  the  coordinate  system  fixed  with  respect  to  the  test  antenna 
and  is  radiating  alone  in  the  homogeneous  fluid  (the  test  antenna  is  removed).  The  expression 
(2.11)  gives  the  far-field  pattern  of  the  test  antenna  as  the  Fourier  transform  of  the  output 
of  the  probe  divided  by  the  far-field  pattern  of  the  probe. 


2.2  Electromagnetic  Fields 

Consider  Figure  2.1  showing  the  geometry  for  measuring  the  electric  field  E  radiated  into 
free  space  by  the  test  antenna  which  is  located  in  the  half  space  z  <  zo  and  is  fed  by  a  signal 
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of  amplitude  cq.  When  the  test  antenna  is  the  only  object  in  all  space  it  radiates  the  field 
E  given  by  the  plane-wave  expansion  [3],  (4,  eqs.(2.10),(2.11)] 

an  r+oo 

E{f)  =  ^  /  ^>^0  (2.12) 

Z'K  »/— 00  J^oo 

where  Tq  is  the  transmitting  spectrum  of  the  test  antenna  given  by 

-  /+00  r+oo  _ 

To{kj;,ky)  =  - -  /  /  z  >  Zq  (2.13) 

ZTTCLq  J—oo  oo 

with  f  =  XX  +  yy  +  zz  and  fc  •  To  =  0.  The  far  field  of  the  test  antenna  is  given  by  the  well 
known  expression  [3],  [4,  eq.(2.50)] 

E{f )  ^ —ioiik  cos  6 - ro(A:cos  (^sin  fcsin  ^sin  0),  z  >  zq  (2.14) 

r 

and  the  far-held  pattern 

^o{6,(i>)  = —ik  cos  6To{k  cos  (l> sin  9,  k  sin  (i>sm  6),  0  <  9  <r  /2.  (2.15) 

for  the  test  antenna  is  defined  such  that  the  far  electric  field  is  E  ao:Fo{9,<f>)e''‘^Jr. 

In  the  following  we  determine  the  spectrum  To  by  measuring  the  electric  field  on  the 
plane  z  =  zq  using  a  probe  as  shown  in  Figure  2.1.  As  for  the  acoustic  field,  the  probe  is 
located  by  its  reference  point  fo  as  measured  in  the  (x,y,z)  coordinate  system  of  the  test 
antenna  (see  Figure  2.1).  Again,  we  assume  that  the  distance  between  the  probe  and  the 
test  antenna  is  always  large  enough  for  multiple  interactions  between  the  test  antenna  and 
the  probe  to  be  negligible.  That  is,  the  reflection  in  the  test  antenna,  of  the  field  radiated 
by  the  induced  currents  on  the  probe,  does  not  affect  the  output  of  the  probe. 

The  scalar  output  of  the  probe  is  denoted  by  bp  and  it  proves  convenient  to  define  the 
vector  output  of  the  probe  as  a  vector  sum  of  the  scalar  output  for  two  mutually  perpendicular 
orientations  [7],  [8]: 

bp  =  6*x  -1-  b^y.  (2.16) 

Here  and  are  the  scalar  outputs  of  the  probe  when  some  reference  line  attached  to  the 
probe  is  parallel  to  the  x  and  y  axis,  respectively.  For  example,  if  the  probe  is  an  open-ended 
rectangular  waveguide  the  reference  line  could  be  chosen  such  that  it  is  parallel  to  the  longest 
cross-sectional  side  of  the  waveguide. 

The  dyadic  receiving  characteristic  of  the  probe  Rp{kx,ky)  is  defined  such  that  the  vec¬ 
tor  output  of  the  probe,  when  it  is  located  at  fo  =  0,  due  to  the  incident  plane  wave 
(27r)"^ee‘^*'*®'*'*»*''*’'’'*^  with  e  •  fc  =  0,  is  given  by  bp  =  Rp{kx,  ky)  ■  e.  Consequently,  the  vector 
output  of  the  probe,  when  it_is  located  at  f©  =  xqX  -t-  yoy  -f  Zqz,  due  to  the  same  incident 
plane  wave,  is  given  by  bp  =  Rp{kx,ky)  •  The  plane-wave  expansion  (2.12) 

for  the  electric  field  radiated  by  the  test  antenna  shows  that  radiated  plane- wave  component 
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propagating  in  the  direction  given  by  {kx,ky)  is  ao(27r)  ^To{kx,  ky)e'^'‘^^'^’‘^^'*"^^'>dkj.dky  and 
thus  the  vector  output  of  the  probe  due  to  this  plane-wave  component  is 

dbyifo)  =  ao%,{kr,ky)  •  foikr,  ky)e'^'‘^^°-^'‘y^°+^^'>dk,dky.  (2.17) 

The  total  vector  output  of  the  probe  is  thus  seen  from  (2.12)  to  be  the  integral  over  the 
entire  kx  —  ky  plane  of  the  expression  in  (2.17),  that  is 

/+00  f  +  OO  —  _ 

/  Rj,{kx,  ky)  ■  fo{kx,  (2.18) 

•oo  J^oo 

which  is  the  final  expression  or  plane-wave  “transmission  formula”  for  the  vector  out 
the  probe  located  at  tq.  Taking  the  inverse  Fourier  transform  of  this  transmission  foi,  la 
shows  that 

Rp{kx,  ky)  ■  fo{kx,  ky)  =  — —  /  /  6^(fo)e-‘(*=*"“+''vvo)rfxodi/o  (2.19) 

(Zir}‘aoJ-oo  J-oo 

from  which  To  will  be  determined.  Since  the  vector  output  of  the  probe  bp  is  normal  to  the 
z  aixis,  (2.19)  gives  only  two  equations  for  determining  the  three  rectangular  components  of 
the  spectrum  Tq.  However,  Tq  satisfies  the  well-known  orthogonality  relation  [4,  eq.(2.12)] 
TQ(kx,  ky)-k  =  0  so  only  two  equations  are  necessary  for  determining  To-  Because  the  vector 
output  of  the  probe  (2.16)  is  perpendicular  to  z  the  receiving  characteristic  can  be  written 
as  _ 

Rp{kx,  ky)  =  iR;{kx,  ky)  -f  0l{kx,  ky)  (2.20) 

with  R^-k  ^0  and  R^-k  =^0.  To  solve  (2.19)  for  the  spectrum  Tq  of  the  electric  field,  define 
the  inverse  of  Rp{kx,  ky)  by  [8] 

ky)  ■  %,{kx,  ky)  =  9{k^,  ky)9{kx,  ky)  +  ^{kx,  ky)^{kx,  ky)  (2.21) 


where  9  and  6  are  unit  vectors  orthogonal  to  k  satisfying  the  relations  9x^  =  k,  kx9  =  ^, 

AAAA  A_  AA 

and  <i>  X  k  =  9  with  k  defined  such  that  kk  =  k  =  kxX  +  kyy  +  jz.  Specifically,  9  and  ^  can 
be  expressed  in  terms  of  kx  and  ky  as 


*’  hJhl  +kl  *  * 


(2.22) 


<f>{kxj  ky)  - 


fkl  +  kl 


fkl  +  kl' 


(2.23) 


If  kx  =  k  cos  (j)  sin  9  and  ky  =  fcsin^sin^,  then  7  =  kcos9,  and  the  unit  vectors  9{kxiky) 

A 

and  4>{kxiky)  are  simply  the  usual  angular  unit  vectors  of  the  spherical  coordinate  system 
{k,9,<i>)  [9,  app.l,  eqs. (115-116)]. 


10 


Writing  R^,  as  ax  +  and  assuming  that  k  ■  x  ^pj  ^  0  one  finds  that  the  inverse 
receiving  characteristic  is  given  by 


=-i  _  fc  X  [R^x  -  R^y 

~  ib  •  [^  X  R^] 


(2.24) 


which  expresses  the  inverse  receiving  characteristic  of  the  probe  in  terms  of  the  vectors 
R^  =  X  ■  Rp  and  R^  =  y  ■  Rp.  Multiplying  (2.19)  with  R^,  and  using  the  definition  of  the 
inverse  receiving  characteristic  (2.21)  show  that  the  spectrum  To  is  given  by 


7'o(fcx,  kp)  =  -nr^Rp  (^r,  ky)-  /  6p(fo)e-‘<"*"‘>+''-'>">)dxodyo.  (2.25) 

(Ziryoo  J-oo  J-oo 

This  is  the  final  probe-corrected  formula  giving  the  spectrum  To  for  the  electric  field  radiated 
into  free  space  by  the  test  antenna  as  a  product  of  the  inverse  receiving  characteristic  of  the 
probe  and  the  Fourier  transform  of  the  vector  output  of  the  probe. 


2.2.1  Formulas  for  Reciprocal  Probes 

Now  consider  the  special  case  of  a  reciprocal  probe  that  is  fed  by  a  waveguide  supporting 
just  one  propagating  mode.  Then  it  follows  from  [3,  ch.2,  eq.(1.6-20a)]  that  the  transmitting 
spectrum  of  the  probe  Tp  is  related  to  the  vector  receiving  characteristic  Rp  by  the  reciprocity 
relation 

Rp{kx,ky)  =  J^-^j^Tp{—kx,—ky).  (2.26) 

Here  c,  p,  and  rjo  are  the  permittivity  of  free  space,  the  permeability  of  free  space,  and  the 
characteristic  admittance  for  the  propagating  mode  in  the  waveguide  feed.  In  Kerns’  work 
[3],  and  So2  correspond  to  our  Tp  and  Rp,  and  sio  corresponds  to  To-  The  vector  receiving 
characteristic  Rp{kx,ky)  of  the  probe  is  defined  according  to  [3,  ch.2,  eqs.(  1.3-8), (1.6- 14a)] 
such  that  the  scalar  output  of  the  probe  (when  it  is  located  at  fo  =  0),  due  to  the  incident 
plane  wave  (27r)"*ec’^*'*^*»‘''*"’'*)  with  e  -  ^  =  0,  is  given  by  bp  =  Rp{kx,  ky)  e.  Letting  T®  and 
T^  be  the  transmitting  spectra  of  the  probe  when  it  is  located  at  Po  =  6  and  its  reference 
line  is  parallel  to  the  x  and  y  axis,  respectively,  one  sees  from  (2.20)  that  R^  and  R^  are  the 
corresponding  vector  receiving  characteristics  of  the  probe.  Applying  the  reciprocity  relation 
(2.26)  to  these  x  and  y  receiving-transmitting  characteristics  gives 

=  =  (2.27) 

One  can  also  see  that  R^  and  R^  are  vector  receiving  characteristics  corresponding  to  in 
Kerns’  work  [3]  by  comparing  the  x  and  y  components  of  the  vector  transmission  formula 
(2.18)  with  Kerns’  scalar  transmission  formula  [3,  ch.2,  eq.(3.1-l),  and  ch.3,  eq.(l-12)]. 
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Equations  (2.27)  allow  the  dyadic  receiving  characteristic  Rj,  in  (2.20)  to  be  written  in 
terms  of  the  transmitting  spectrum  of  the  probe  in  two  different  orientations.  One  may  also 
insert  the  values  of  f^{kx,ky)  and  ^{kx,ky)  from  (2.27)  into  the  expression  (2.24)  for  the 

inverse  of  Rp  and  then  use  (2.25)  to  calculate  the  transmitting  spectrum  of  the  test  antenna. 

Let  us  finally  express  the  far-field  pattern  (2.15)  of  the  test  antenna  in  terms  of  the  vector 
output  of  the  probe  and  the  far-field  patterns,  and  of  the  probe  in  the 

two  orientations  when  it  is  located  at  fo  =  6.  To  do  this,  we  use  the  far-field  relation  (2.14) 
and  the  expression  (2.25)  for  the  transmitting  spectrum  to  show  that  the  far-field  pattern 
for  the  test  antenna  is  given  by 

-  ,  ik  cos  0  =-l  f+oo  r+oo  _  _ 

R  {kcos4>sm6,ksm<f>sin6)  ■  /  6p(ro)e~‘*'’’’’®dxodyo-  (2.28) 

J-oo  J-oa 

Since  the  far-field  pattern  of  the  probe  is  to  be  evaluated  in  the  hemisphere  with  z  <  zq  \t  is 
found  that  [3,  ch.2,  eq.(1.2-16b)],  [4,  sec.2.2.1]  =  ifccos^T* (A:cos  (^sin0,  Aisin  <^sin^) 

for  7r/2  <  6  <ir,  and  the  reciprocity  relation  (2.27)  then  shows  that 


Rp{k  cos  ^  sin  ^ ,  fc  sin  <f>  sin  B) 


^^{■K  -  6l,7r  -I-  <^) 


(2.29) 


and  similarly  for  R^.  Inserting  these  results  into  the  expression  for  the  inverse  receiving 
characteristic  (2.24)  yields 


itfokf  X 

-9,ir  +  <l>)x-  -9,n  + 

•f>)y 

^®(7r  -  ^,7r  +  ^)  X  ^{ir  -  9,ir  +  ^)j 

(2.30) 

which  is  the  final  expression  giving  the  inverse  of  the  dyadic  receiving  characteristic  of  the 
probe  in  terms  of  the  far-field  patterns  of  the  probe  in  two  orientations.  Inserting  (2.30)  into 
(2.28)  gives  the  far-field  pattern  of  the  test  antenna  in  terms  of  the  far-field  patterns  of  the 
probe  and  the  Fourier  transform  of  the  vector  output  of  the  probe 


fE  f  X  -  g,ir  -b  <i>)x  -  -B,ir-\-  <l>)y\ 


(2.31) 

J-^QO  J^oo 

If  the  receiving  probe  is  a  reciprocal  elementary  electric  dipole,  the  formula  (2.31)  for 
the  far  field,  or  more  generally,  the  formula  (2.25)  for  the  plane-wave  spectrum  of  the  test 
antenna  should  reduce  to  the  non-probe-corrected  formulas  of  [4].  To  prove  this,  begin  with 
the  expression  for  the  receiving  characteristic  and  output  of  a  reciprocal  electric  dipole  probe 
(3,  p.131] 

R;  =  X  (it  X  x)  (2.32) 

^  AirtooTio  ''  ' 
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^  X  »)  (2.33) 

^  =  +  =  (2.34) 

where  p  is  the  magnitude  of  the  electric  dipole  moment  and  Exy  =  ExX  +  EyP  is  the  transverse 
electric  field  of  the  test  antenna  in  the  scan  plane.  Substituting  and  from  (2.32)  and 

(2,33)  into  (2.24),  and  simplifying  the  resulting  dyadic  expression  shows  that  *  for  the 
dipole  probe  can  be  written  as 


2^r  - 

Rp  = - k  X  (yx  -  xy) .  (2.35) 

7 

With  bp  and  R^  from  (2.34)  and  (2.35)  inserted  into  (2.25),  the  spectrum  of  the  test  antenna 
takes  the  form 


_  _  -  /+00  r+oo 

Tikx,  ky)  =  aoToikx,  ky)  =  — - - k  X  /  i  X  .e,„(fo)e-'(***®+''»"®)dxodyo  (2.36) 

ZTT'Y  oo  oo 


or 


fikx,ky) 


0— t-yzo  f+oo  f+oo 


2ir 


n 

J—OO  •/  — OO 


Exy{ro)  ~  ^zk  ■  Exy{fo) 


(2.37) 


Equation  (2.37)  determines  the  plane- wave  spectrum  of  the  test  antenna  in  terms  of  the 
transverse  electric  field  of  the  test  antenna  on  the  scan  plane.  It  is  identical  to  the  result 
obtained  for  the  spectrum  of  the  test  antenna  from  the  non-probe-corrected  formulas  [4, 
eqs.(2.11)-(2.12)]. 

In  the  next  chapter,  the  frequency-domain  formulas  of  this  chapter  for  probe-corrected 
near-field  measurements  will  be  translated  into  the  time  domain. 
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Chapter  3 

Time-Domain  Formulas 


The  probe- corrected  planar  near-field  time-domain  formulas  are  derived  in  this  chapter  for 
both  acoustic  and  electromagnetic  fields.  Section  3.1  derives  the  formulas  for  the  acoustic 
field  first,  by  taking  the  inverse  Fourier  transform  of  the  corresponding  frequency-domain 
formulas  of  Section  2.1,  and  then  directly  in  the  time  domain.  Section  3.2  derives  the 
corresponding  formulas  for  the  electromagnetic  field.  Finally,  Section  3.3  uses  the  probe- 
corrected  formulas  to  numerically  calculate  the  far-field  pattern  of  an  acoustic  point  source 
from  near-field  data  obtained  by  a  nonideal  probe. 


3.1  Acoustic  Fields 

Without  loss  of  generality  we  assume  that  the  scan  plane  is  at  z  =  0  because  it  simplifies 
the  time-domain  analysis  of  this  section  significantly.  Thus,  the  test  antenna  is  located  in 
the  half  space  z  <  0  and  radiates  the  time-domain  acoustic  field  $(r,  t)  satisfying  the  scalar 
wave  equation 

=  j:>0  (3.1) 

where  c  is  the  signal  speed.  The  acoustic  field  $(f,  t)  can  be  found  from  its  Fourier  transform 
$„(f)  given  in  Section  2.1,  by  the  Fourier  inversion  formula 

/+00  .  ^ 

^u,(r)e-’"‘dw.  (3.2) 

'OO 

In  this  chapter  all  frequency-domain  fields  are  labeled  with  subscript  lj.  The  time-domain 
probe-corrected  formulas  are  derived  in  Subsection  3.1.1  using  the  Fourier  inversion  formula 
(3.2)  and  the  frequency-domain  probe-corrected  formulas  of  Section  2.1.  Subsection  3.1.3 
rederives  these  formulas  by  working  directly  in  the  time  domain  with  the  Radon  transform 
and  a  time-domain  receiving  characteristic  for  the  probe.  It  is  shown  that  the  probe-corrected 
formulas  simplify  significantly  when  the  probe  output  due  to  an  incoming  plane  wave  is 
proportional  to  the  time  derivative  of  the  field  of  that  plane  wave. 


Before  the  time-domain  probe-corrected  formulas  are  derived  we  summarize  the  time- 
domain  field  representation  from  [4,  sec.3.2.2)  on  which  all  the  time-domain  formulas  of  this 
report  are  based.  We  start  by  defining  the  two  time-domain  functions 


ro{Cri,t)=  I  To^{uj^,u;r})e 

•/  — OO 


/+O0  .  1  z+cx) 

aou,To^{u}^,  <^  =  /  ao{t-  t')To{(,  t],  i')di'  (3.4) 

where  Tow  is  the  plane-wave  spectrum  (2.2)  of  the  test  antenna,  and  (^,  t])  are  new  spectral 
variables,  and  Oow  is  the  spectrum  of  the  input  signal  such  that  ao(<)  =  aQ^e~"^^du)  is 

the  time  dependence  of  the  input  signal  to  the  test  antenna.  Use  the  definition  of  T  along 
with  the  expression  (2.2)  for  the  frequency-domain  spectrum  to  get  [4,  eq.(3.51)] 

2  f+oo  r+oo 

»/.  0  =  TT  /  /  ^(^0,  t  +  (xo  +  T}yo)dxodyo,  ro  =  Xqx  yoy  (3.5) 

which  shows  that  T  is  simply  the  Radon  transform  of  the  time-domain  field  $(fo,<)  on  the 
plane  z  =  0  when  the  test  antenna  is  fed  by  the  input  signal  ao(t).  The  results  of  [4]  show 
that  the  time-domain  field  in  the  half  space  z  >  0  is  given  by  [4,  eq.(3.56)] 

wn^r#» 

.  r  e  +  ri^<c-^ 

^  -  \  iWF+^r^u  e+r}^>  c-^  •  ^  ^  ^ 


^  -  \  iWF+^r^u  e+r}^>  c-^  •  ^  ^  ^ 

Furthermore,  the  time-domain  far  field  is  given  by  the  simple  expression  [4,  eq.(3.59)l 

$(f,  <)  ~ ———T(c“*  cos  ^  sin  ^,c~’ sin  ^  sin  —  r/c),  0<9<ir/2.  (3.8) 

cr  at 

These  expressions  show  that  the  field  radiated  by  the  test  antenna  into  the  half  space  z  >  0 
is  completely  determined  from  the  Radon  transform  T{(,r],t)  of  the  field  on  the  plane  z  =  0. 
Writing  the  time-domain  far-field  as  $(f ,  t)  ~  ^{0,  <f>,t  —  rf  c)/r,  where  T  is  the  time-domain 
far-field  pattern,  we  find  from  (3.8)  that 

A  Q 

.F(0,  t)  =  — — ■5-T(c~*  cos<^sin61,c~’ sin<^sin^,t),  0  <  0  <  7r/2.  (3.9) 

c 

The  formulas  (3.5)-(3.9)  constitute  the  time-domain  field  representation  that  gives  the  field 
in  the  half  space  z  >  0  in  terms  of  the  Radon  transform  of  the  field  on  the  plane  z  =  0. 

Having  summarized  the  time-domain  field  representations,  we  will  now  derive  expressions 
that  give  the  Radon  transform  T  or  the  far-field  pattern  T  in  terms  of  the  output  and 
receiving  characteristic  of  the  probe. 
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3.1.1  Formulas  Obtained  from  the  Fourier  Transform 

We  begin  by  deriving  the  time-domain  probe-corrected  formulas  using  the  Fourier  transform 
along  with  the  frequency-domain  formulas  of  Section  2.1.  The  frequency-domain  probe- 
corrected  formula  (2.8)  shows  that  with  fo  =  xqx  -f  yoy 

(3.10) 

where  Rj^  and  bp^  are  the  frequency-domain  receiving  characteristic  and  output,  respectively, 
of  the  probe.  Introduce  the  inverse  time-domain  receiving  characteristic  of  the  probe  defined 

by’ 


1  /+00 
•/— oo 


O— tu/t 


oo 


du>. 


Use  the  convolution  rule 


/*”  /^“  /((')»((  -  t')<U' 

J^oo  Zir  •/— oo 


(3.11) 


(3.12) 


in  taking  the  inverse  Fourier  transform  of  (3.10)  to  get  the  following  three  alternative  formulas 
for  the  Radon  transform  T 

/  Rp  ((^Vii)  /  bpifo,t  -  t  -i- ^xo  +  vyo)dxodyodt'  (3.13) 

\Z‘K)'^J~oo  J~oo  J—oo 

1  f+OO  f+OO  f+OO 

/  /  /  Rp  {(,r!,t-t+(xo  +  riyo)bp{foj')dxodyodt'  (3.14) 

./-oo  J—oo  J—oo 

and 

1  r+oo  r+oo  f+OO 

7’(^,»7,t)  =  7;r-^  /  /  /  Rp  i^,V,*  -  ^  +  (xo  + vyo)bpifo,i')dt'dxodyo  (3.15) 

yZir)  J  —  oo  J  —  oo  J  —  oo 

where  fo  =  XoX  +  yoy  and  6p(fo,  t')  is  the  time-domain  output  of  the  probe  when  it  is  located 
at  ffl.  The  formulas  (3.13),  (3.14),  and  (3.15)  are  the  general  time-domain  probe-corrected 
formulas  giving  the  Radon  transform  of  the  field  in  terms  of  the  output  and  receiving  char¬ 
acteristic  of  the  probe. 

Assume  that  the  test  antenna  is  turned  on  at  f  =  to  so  that  all  fields  are  zero  for 
t  <  to-  Then  the  output  of  the  probe  6p(fo,<)  is  zero  for  fo  ^  p{t)  where  p(t)  is  the  finite 
region  of  the  scan  plane  on  which  the  field  is  nonzero  at  time  t.  (For  late  times,  p(t)  is 
approximately  a  circular  disk  with  radius  c(t  —  to)-)  Without  making  any  assumption  about 
the  time  dependence  of  the  inverse  receiving  characteristic  the  probe-corrected  formulas 
all  involve  at  least  one  integral  over  an  infinite  region.  For  example,  the  (lojJ/o)  integral  in 

*Both  the  time-domain  receiving  characteristic  Rp(^,r),t)  =  Rpu(u>(,iJf))e~"^*du  and  its  inverse 
defined  by  (3.11)  can  be  generalized  functions.  This  will  be  seen  from  Section  3.1.3  where  these  functions 
are  defined  directly  in  the  time  domain. 


(3.14)  becomes  an  integral  over  the  finite  region  /j(t'),  but  the  t'  integration  goes  from  t'  =  to 
to  t'  =  +00  and  thus  involves  an  integration  over  an  infinite  region.  Similar  statements  hold 
for  the  probe- corrected  formulas  (3.13)  and  (3.15).  However,  as  we  shall  see  later  in  this 
section,  only  integrals  over  finite  regions  need  be  calculated  to  get  the  time-domain  far- field 
pattern  when  a  special  type  of  probe  is  used. 

Equation  (3.15)  may  be  the  most  attractive  for  measurement  purposes  since  the  inner 
time  integral  can  be  calculated  at  each  scan  point  fo  for  the  desired  values  of  Thus 

one  would  not  have  to  store  the  values  of  the  output  of  the  probe  6p(fo,  <')  over  time  to 
calculate  the  Radon  transform  of  the  field  radiated  by  the  test  antenna. 

Having  determined  the  Radon  transform  T{(,T],t)  one  can  then  use  (3.6)  to  calculate 
the  time-domain  field  everywhere  in  the  half  space  z  >  0  or  use  (3.9)  to  calculate  the  time- 
domain  far-field  pattern.  Note  also  that  instead  of  the  real  time-domain  Radon  transform 
formuleis  (3.5)-(3.8)  one  could  have  used  the  complex  formulas  of  (4,  sec. 3. 2. 3]  involving  the 
analytic  fields.  The  expressions  for  the  field  in  the  half  space  s  >  0  would  then  have  been 
given  by  the  simpler  expression  [4,  eq.(3.78)],  which  however  involves  the  analytic  fields  that 
cannot  be  measured  directly.  Thus,  for  measurement  purposes,  the  real  Radon  transform 
formulas  used  in  this  report  seem  advantageous  compared  to  the  complex  Radon  transform 
formulas.  This  was  discussed  more  fully  in  the  previous  report  [4]. 


3.1.2  Formulas  for  Time-Derivative  Probes 


The  general  time-domain  probe- corrected  formulas  (3.13)-(3.15)  require  the  calculation  of  a 
triple  integral  involving  an  infinite  time  integral.  We  will  now  show  that  for  a  special  type 
of  probe,  which  has  been  used  for  time-domain  electromagnetic  measurements  [10],  [11], 
the  infinite  time  integral  can  be  evaluated  analytically.  The  outputs  of  these  probes  due  to 
incoming  plane  waves  are  proportional  to  the  time  derivative  of  the  field  in  the  incoming 
wave,  that  is,  the  frequency-domain  receiving  characteristic  is  proportional  to  u;  and  can  be 
written  as 

(3.16) 

where  the  angularly  dependent  function  Qp{^,  fj)  is  independent  of  the  frequency.  With  the 
receiving  characteristic  (3.16)  inserted  into  the  probe-corrected  formula  (3.10),  the  inverse 
Fourier  transform  of  (3.10)  gives 

d  1  /■+°°  /■+°° 

=  2TfQ  r))J  oo  J-oo  +  + 

which  expresses  the  time  derivative  of  the  Radon  transform  of  the  near  field  as  the  Radon 
transform  of  the  output  of  the  probe  divided  by  the  angularly  dependent  function  Qp.  The 
far-field  expression  (3.9)  shows  that  the  time-domain  far-field  pattern  is  given  by 


cosg  bp{ro,t  +  r  ■  rolc)dxodyQ 

27rcQp(c“*  cos  ^ sin  6,  c~^  sin  4>  sin  0) 


(3.18) 


Thus,  when  the  probe  has  the  simple  receiving  characteristic  (3.16),  the  far  field  of  the  test 
antenna  is  simply  the  Radon  transform  of  the  output  of  the  probe  divided  by  the  angular 
dependence  of  the  probe  receiving  characteristic. 

Let  us  again  assume  that  the  test  antenna  is  turned  on  at  <  =  <o  so  that  the  output 
of  the  probe  bp{fo,t)  is  zero  when  fo  is  outside  the  region  p{t).  Then  the  integrand  of 
(3.17)  is  zero  when  fo  is  outside  the  region  p{t  +  ^Xq  +  VVo)-  For  late  times,  when  the  finite 
region  p{t)  is  approximately  a  circle  of  radius  c{t  —  to),  this  translates  into:  6p(fo,t)  =  0  for 
Ifol  >  c{t  -  to)  +  c^xo  +  CTjyo.  When  +  t)^  >  c~^  this  condition  is  not  satisfied  no  matter 
how  large  |fo|  may  be,  so  in  this  situation  the  region  of  integration  in  (3.17)  is  infinite.  When 
this  condition  is  satisfied  for  all  |fo|  larger  than  some  finite  number  and  the 
region  of  integration  in  (3.17)  (and  (3.18))  is  finite.  Consequently,  when  the  time-domain 
far-field  pattern  is  calculated  with  a  probe  whose  receiving  characteristic  satisfies  (3.16),  the 
integral  is  over  a  finite  region.  Thus,  the  far-field  pattern  can  be  calculated  at  early  times 
from  measured  data  taken  on  a  finite  scan  plane  without  introducing  any  error  caused  by  the 
finite  scan  plane.  This  property  of  the  special  time-domain  probe- corrected  formula  (3.18) 
is  not  possessed  by  any  of  the  frequency-domain  probe-corrected  formulas  and,  as  discussed 
in  [4,  ch.4],  makes  it  possible  to  use  planar  time-domain  scanning  for  broad  beam  antennas. 


3.1.3  Formulas  Obtained  Directly  in  the  Time  Domain 

This  subsection  derives  the  probe- corrected  acoustic  formulas  directly  in  the  time  domain 
using  the  Radon  transform  and  a  time-domain  receiving  characteristic  for  the  probe.  To 
define  a  time-domain  receiving  characteristic  of  the  probe,  one  has  to  specify  the  time-domain 
basis  fields  that  will  be  used  to  expand  the  field  radiated  by  the  test  antenna.  We  choose  the 
time-domain  basis  fields  to  be  simply  the  inverse  Fourier  transform  of  the  propagating  and 
decaying  plane  waves  used  to  expand  the  frequency- domain  fields.  To  facilitate  the  time- 
domain  analysis,  the  spectral  variables  {kj„ky)  are  replaced  by  (^,t/)  defined  by  {kx,ky)  = 
{u)^,u}Ti)  for  all  real  u  and  we  only  consider  fields  in  the  half  space  z  >  0.  Furthermore,  we 
write  the  frequency-domain  plane  wave  (2T)"’e’^***''’*'» as 


_1_  r  u;  >  0 


(3.19) 


where  — oo  <  ^  <  -l-oo,  — cx3  <  t]  <  -hoc,  z  >  0,  and  *  denotes  complex  conjugation. 
Taking  the  inverse  Fourier  transform  (3.2)  of  equation  (3.19)  shows  that  the  time-domain 
basis  fields  are 


{ 


6{t  -  (x  -  rjy  -  (z), 
_ _ 


e+r,-‘>c--‘ 


where  z  >  0  and  the  identity 


2|CN 


z  >  0 


(3.20) 


(3.21) 
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Figure  3.1:  Propagating  time-domain  plane  wave. 


has  been  used. 

Since  these  time-domain  basis  fields  are  obtained  from  frequency-domain  plane  waves,  we 
will  refer  to  them  as  time-domain  plane  waves.  As  will  be  seen  below,  the  time-domain  plane 
waves  $®(f,  Tj)  constitute  a  complete  set  of  functions  for  expanding  the  time-domain  field 
radiated  by  the  test  antenna. 

Let  us  study  the  physical  properties  of  these  time-domain  plane  waves.  The  propagating 
plane  wave  <  c~^  is  a  delta  function  which  is  nonzero  only  on  the 

plane  <  =  -f  -f  (z  and  is  propagating  in  the  direction  given  by  the  unit  vector  = 

c^x  +  CTjy  +  c(z  as  shown  in  Figure  3.1.  The  decaying  plane  wave  ^^  +  7^  > 

is  more  complicated.  On  the  plane  ^zd-yy  =  to»  which  is  parallel  to  the  z  axis,  the  amplitude 
of  the  decaying  plane  wave  has  a  maximum  at  z  =  |t  —  to|/|CI  decays  as  z~^  at  infinity  as 
shown  in  Figure  3.2.  Furthermore,  from  the  analysis  of  Morse  and  Feshbach  [13,  pp.813-814] 
it  follows  that 


r-i 


^(Cl 


z|<|™o+  z^lCP  -1-  (t  -  -  rjyy 


=  6{t  -^x  -  Tjy) 


(3.22) 


so  that,  in  agreement  with  (3.19),  all  time-domain  plane  waves  are  delta  functions  on  the 
plane  z  =  0. 

Having  defined  and  studied  the  time-domain  plane  waves,  let  us  now  turn  to  the  time- 
domain  receiving  characteristic  of  the  probe.  This  characteristic  is  denoted  by  Rp{^,  y,  t)  and 
is  defined  to  be  the  output  of  the  probe,  when  its  reference  point  is  located  at  tq  =  0,  and  the 
incoming  wave  is  the  plane  wave  $®(f,t,^,y)  given  in  (3.20).  From  the  expression  (3.20)  for 
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Figure  3.2:  Decaying  time-domain  plane  wave  on  the  plane  -1-  r/y  =  to- 

the  plane  wave  it  is  seen  that  the  output  of  the  probe,  when  it  is  located  at  fo  =  xqx  -f  yoy, 
is  given  by  -  ^lo  -  Wo)  as  illustrated  in  Figure  3.3. 

Now  that  both  the  time-domain  plane  waves  and  the  time-domain  receiving  characteristic 
of  the  probe  have  been  defined,  v.  can  calculate  the  output  of  the  probe  due  to  the  field 
radiated  by  the  test  antenna.  To  do  this  note  that  the  expression  (3.6)  for  the  field  radiated 
by  the  test  antenna  can  be  written  in  terms  of  the  time-domain  plane  waves  as  the  time- 
convolution  integral 

1  r+’x  f+oo  r+oo  fp 

/-CO  Loo  ^  z>0.  (3.23) 

This  expression  enables  us  to  derive  the  expression  for  the  output  of  the  probe  6p(fo,  t)  when 
it  is  located  at  fo  and  the  incident  field  is  that  of  the  test  antenna.  We  assume  that  the 
distance  between  the  probe  and  the  test  antenna  is  large  enough  that  multiple  interactions 
between  the  test  antenna  and  the  probe  can  be  neglected.  That  is,  the  reflection  in  the 
test  antenna,  of  the  field  radiated  by  the  induced  currents  on  the  probe,  does  not  affect  the 
output  of  the  probe. 

The  output  of  the  probe  from  the  component  — (27r)~'$®(f,  t  — t',  y)^T(^,  y,  t')dt'd^d‘q 
of  the  incident  field  is,  according  to  the  definition  of  the  probe  receiving  characteristic,  given 

2 

d6p(fo,0  =  ~  ~  ^xo  -  T]yo)-^T{(,ri,t')dt'd^drj  (3.24) 
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Figure  3.3:  Output  of  probe  due  to  the  incoming  plane  wave  The  reference 

point  of  the  probe  is  located  at  fo  =  Xqx  +  yoy. 

« 

and  the  total  output  of  the  probe  is  thus  seen  from  (3.23)  to  be 

J  f+oo  f+oo  r+oo  ffi 

bp{fo,t)  =  (3.25) 

To  determine  the  Radon  transform  from  this  equation,  define  the  inverse  time- 

domain  receiving  characteristic  for  the  probe  in  accordance  with  (3.11)  by 

/+00 

-  t')dt'  =  47r26(t),  (3.26) 

'OO 

multiply  the  expression  (3.25)  for  the  total  output  of  the  probe  by  R~^(^,i],t"  —  t),  and 
integrate  it  with  respect  to  t  to  get 

/+00  r+oo  r+oo 

^  bp{fo,t)R;\(,ri,t"  -  t)dt  =  -2x  J_^  -  ^xo  -  vyo)d(drj.  (3.27) 

The  Radon  transform  pairs  (3.5)  and  (3.6)  (see  also  [12,  p.lll]) 

1  /+°°  r+oo 

=  y  (xq  +  T]yo)dxodyo  (3.28) 

$(fo,t)  =  y^^  y^^  -^xo  -  T)yo)d(dT] 
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(3.29) 


with  fo  =  XqX  +  yoy,  along  with  (3.27),  now  show  that  the  Radon  transform  of  the  field 
radiated  by  the  test  antenna  on  the  plane  z  =  0  is  given  by 

1  f+oo  r+oo  f+OO 

=  7:7-^  /  /  /  Rp  {i,r},t-t'  +  (xQ  +  T)yo)bp{fQ,t')dt'dxodyo.  (3.30) 

j-00  J-00  J-00 

This  expression  is  seen  to  agree  with  the  result  (3.15)  found  by  Fourier  transforming  the 
corresponding  frequency-domain  result.  Equation  (3.30)  is  the  time-domain  probe-corrected 
formula  derived  directly  in  the  time  domain  using  the  Radon  transform  and  a  time-domain 
expansion  of  the  field  radiated  by  the  test  antenna.  Interchanging  the  orders  of  integration 
in  (3.30)  produces  the  probe- corrected  formulas  (3.13)  and  (3.14). 


Formulas  for  time>derivative  probes 

Let  us  again  turn  to  the  special  case  where  the  output  of  the  probe  due  to  an  incoming 
plane  wave  is  proportional  to  the  time  derivative  of  the  plane  wave.  In  this  case  the  time- 
domain  receiving  characteristic  takes  the  form 

V,  t)  =  Q,((,  t,(,  n)  =  QA(<  iW)  (3.31) 

which  is  valid  for  all  (^,  J7)  and  is  in  accordance  with  the  frequency-domain  result  (3.16).  The 
last  part  of  (3.31)  is  obtained  by  using  the  expressions  (3.20)  and  (3.22)  to  get  the  values 
of  the  plane  waves  on  the  plane  z  =  0.  The  definition  (3.26)  for  the  time-domain  inverse 
receiving  characteristic  and  the  expression  (3.31)  for  the  receiving  characteristic  show  that 


Qr{(.V) 


(3.32) 


for  all  (^,7).  Taking  the  time  derivative  of  the  probe-corrected  formula  (3.30)  and  inserting 
this  value  for  the  time  derivative  of  R~^,  one  recovers  the  probe-corrected  formula  (3.17) 
derived  from  the  corresponding  frequency-domain  result.  It  has  now  been  demonstrated 
once  more  that  the  probe-corrected  formulas  simplify  considerably  when  the  output  of  the 
probe  due  to  an  incoming  plane  wave  is  proportional  to  the  time  derivative  of  the  plane 
wave. 


3.1.4  Formulas  for  Reciprocal  Probes 

In  the  previous  two  subsections  the  probe-corrected  formulas  were  derived  for  arbitrary 
probes  and  for  probes  whose  output  is  proportional  to  the  time  derivative  of  an  incoming 
plane  wave.  This  subsection  derives  probe-corrected  formulas  for  cases  wh**"e  the  probe  is  a 
reciprocal  electroacoustic  transducer  [6].  It  is  shown  that  there  is  a  simpk  relation  between 
the  time-domain  receiving  characteristic  Rp{(,r},t)  and  the  time-domain  transmitting  char¬ 
acteristic  Tp{^,rf,t)  of  the  probe.  Here  is  simply  the  Radon  transform  (3.5)  of  the 

field  radiated  by  the  probe  into  the  half  space  z  <  0  when  it  is  located  at  fp  =  0  and  is  fed 
by  the  delta-function  input  signal  ap(t)  =  2ir6{t),  as  shown  in  Figure  3.4. 
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Figure  3.4:  Geometry  for  reciprocity  relations. 


The  field  radiated  by  the  probe  into  the  half  space  z  <  0  in  this  situation  is  denoted  by 
and  is  given  by  the  time-convolution  integral  (3.23)  with  T  replaced  by  T^,  and  with 
being  the  time-domain  plane  waves  for  the  half  space  z  <  0  given  by  (3.20)  with  z  replaced 
by  —z.  The  delta-function  far-field  of  the  probe  is  given  by  the  negative  of  (3.8),  that  is 

A  Q 

~  ~~~ ‘HI ^»(^~*  sin  ^,c~*  sin  sin  —  r/c),  7r/2  <  0  <  tt  (3.33) 

^  cr  ot  ^ 

and  the  corresponding  far-field  pattern  is 

^(6, cos  <f>sm0yC~^  sin  <^sin0,t),  ■k/2<0<-k  (3.34) 

^  c  at 

where  is  defined  such  that  the  delta- function  far-field  is 

Assume  now  that  the  probe  is  a  reciprocal  electroacoustic  transducer  [6].  Then  from 
(2.9)  the  following  frequency-domain  reciprocity  relation  exists  for  all  real  u  between  its 
frequency-domain  receiving  and  transmitting  characteristics 


(3.35) 


where  the  fact  that  any  frequency-domain  function  must  satisfy  has  been  used. 

Taking  the  inverse  Fourier  transform  of  (3.35)  and  using  the  fact  that  the  real  and  imaginary 
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parts  of  an  analytic  function  on  the  real  axis  are  related  through  the  Hilbert  transform  [13, 
pp. 370-372]  reveal  the  following  time-domain  reciprocity  relation 


(3.36) 


(3.37) 


where  H  denotes  the  Hilbert  transform  given  by 

/■+«> 

The  bar  on  the  integral  sign  indicates  a  Cauchy  principal  value  integration  (see  [4,  sec.3.2.3]). 
Note  that  the  time-domain  reciprocity  relation  (3.36)  is  simple  and  similar  to  the  correspond¬ 
ing  frequency-domain  relation  (2.9)  as  long  as  +  <  c~^  (propagating  plane  waves);  how¬ 
ever  it  is  more  complicated  and  involves  the  Hilbert  transform  when  +  >  c~^  (decaying 

plane  waves). 

The  time-domain  reciprocity  relation  (3.36)  with  C,  —  — c~’  cos0  now  gives  the  following 
relation  between  the  delta-function  far-held  pattern  (3.34)  and  the  receiving  characteristic 
of  the  probe 

^(0,  <f>,  t)  =  TfQpQ—Rp{c~^  cos(7r  -I-  (j>)  sin(7r  —  6),  c~’  sin(7r  -f-  <f>)  sin(7r  —  0),t),  irl2  <  0  <  ir. 

(3.38) 

Consequently,  for  reciprocal  probes,  the  delta-function  far-field  pattern  is  proportional  to  the 
time  derivative  of  the  time-domain  receiving  characteristic.  Since  the  delta-function  far-held 
pattern  is  equal  to  the  time  derivative  of  the  step-function  far-held  pattern,  it  is  seen  that 
the  step-function  far-held  pattern  of  a  reciprocal  probe  is  proportional  to  its  time-domain 
receiving  characteristic. 

Finally  assume  that  the  probe  receiving  characteristic  satishes  (3.31)  in  addition  to  the 
reciprocity  relation  (3.36),  that  is,  the  output  of  the  probe  from  an  incoming  plane  wave  is 
proportional  to  the  time  derivative  of  the  plane- wave  held.  Then  (3.31)  combines  with  (3.38) 
to  produce 

J^{0, 4>,  i)  =  rioPoQv{c~^  cos(7r  -f  4>)  sin(7r  —  0),  c~^  sin(7r  -|-  ^)  sin(7r  -  0))S"{t),  7r/2  <  0  <  it 

(3.39) 

which  shows  that  the  angular  dependence  of  the  delta-function  far-held  is  that  of  the  receiving 
characteristic  and  the  time  dependence  of  the  far-held  pattern  is  the  second  derivative  of  a 
delta  function.  Consequently,  the  far-held  of  this  particular  reciprocal  probe  is  related  to 
the  second  time  derivative  of  the  signal  feeding  the  probe. 


3.2  Electromagnetic  Fields 

This  section  derives  the  probe-corrected  formulas  for  the  electromagnetic  helds.  The  deriva¬ 
tions  are  very  similar  to  the  ones  for  the  acoustic  held  since  each  of  the  rectangular  com¬ 
ponents  of  the  electric  held  satishes  the  acoustic-held  formulas.  Therefore,  nearly  all  of  the 


formulas  for  the  electromagnetic  fields  can  be  obtained  easily  by  modifying  the  corresponding 
formulas  for  the  acoustic  field. 

First  we  summarize  the  Radon  transform  formulas  for  the  electric  field,  which  allow  us 
to  calculate  the  electric  field  everywhere  in  the  half  space  2  >  0  in  terms  of  the  field  on  the 
plane  2  =  0.  We  assume  that  the  test  antenna  is  fed  by  the  time  signal  ao(0  and  introduce 
the  Radon  transform  of  the  electric  field  on  the  plane  2  =  0  by  [4,  eq.(3.61)] 

I  r+00  r+00  _  _ 

=  —  /  /  E(ro,t  +  (xo  +  Tjyo)dxodyo,  tq  =  XqX  +  yoy.  (3.40) 

Z'TT  00  00 

Because  V  •  jE  =  0  the  Radon  transform  of  the  electric  field  satisfies  [4,  eqs.(3.64),(3.67)] 

+  +  =  +  (3.41) 

ru,  V,t)Ux-\-  rjy)  +  \C\nm,  r,,t)-z  =  0,  (3.42) 

where  Ti  is  the  Hilbert  transform  given  in  (3.37).  Since  V  x  £  =  it  follows  that  the 

Radon  transform  for  the  magnetic  field  is  [4,  eqs.(3.69),(3.70)] 

ThK. >).<)  =  +  +  {’  +  >?“<<:■*  (3.43) 

?«(«,  <?,«)  =  +  to)  X  ?(<.>).<)  +  — 5  +  (3.44) 

fi  n 

and  thus  the  Radon  transform  of  the  magnetic  field  is  determined  from  the  Radon  transform 
of  the  electric  field. 

The  electric  field  in  the  half  space  2  >  0  can  be  calculated  by  the  formula  [4,  eq.(3.62)] 

I  f+00  f+00  f+00  - 

“  “s /-„  /-o,  *  ^ 

where  is  the  time-domain  plane  wave  (3.20).  The  time-domain  far  field  is  given  by  the 
simple  expression  [4,  eq.(3.71)] 

E{f,t) cos<f)sin6,c~^  s\nd>s\n0,t  —  r/c),  O<0<n/2.  (3.46) 

cr  at 

Defining  the  time-domain  far-field  pattern  for  the  test  antenna  such  that  the  far 

electric  field  is  E{f,t)  ~  —  rjc)lr  one  finds  from  (3.46)  that 

^{0,<l>,t)  =  cos<l>sm0,c~^  s'm<i>sin0,t),  O<0<ir/2.  (3.47) 

c  di 

Having  summarized  all  the  Radon  transform  formulas  for  the  electric  field,  we  now  turn  to 
the  problem  of  deriving  the  time-domain  probe-corrected  formulas.  As  in  Section  2.2  for  the 
frequency-domain  electromagnetic  field,  the  vector  output  of  the  probe  is  defined  by  (2.16) 
as  a  vector  sum  of  the  scalar  output  in  two  orientations.  The  time-domain  dyadic  receiving 


characteristic  Rp{(,ri,t)  is  defined  such  that  when  the  incident  field  is 

with  e(^,T/)  satisfying  (3.41)  and  (3.42),  the  vector  output  of  the  probe  is  given  by 

6p(fo,t)  = /2p(^,i7,t  -  ^lo  -  i?l/o)  •  e(^,»/),  ro  =  iox  +  i/oy  (3.48) 


when  the  probe  is  located  at  tq. 

The  time-domain  dyadic  receiving  characteristic  is  simply  the  inverse  Fourier  transform 
of  the  corresponding  frequency-domain  dyadic  and  because  the  vector  output  of  the  probe 
has  only  x  and  y  components,  it  can  be  written  as 

RpiiyVJ)  =  +  yR^{^,ri,t).  (3.49) 

From  the  expansion  (3.45)  of  the  field  radiated  by  the  test  antenna  it  is  seen  that  when  the 
probe  is  illuminated  by  this  field  and  is  located  at  fo,  its  total  output  is  given  by 

1  r+o°  /+O0  f+oo—  Q'i  _ 

6p(ro,0  = y_^  y_^  Rp{(yVJ-i'-^xo-T]yo)  —T{^,T),t')dt'd^dTi  (3.50) 


which  is  the  electromagnetic  analog  to  the  acoustic  formula  (3.25).  The  time-domain  inverse 

receiving  characteristic  R^  is  simply  the  inverse  Fourier  transform  of  the  frequency- 

domain  inverse  receiving  characteristic  given  in  (2.24) 


ktv)  X  -  R^{uj^,u^T))y] 

kU,r})  •  X  fi^(w^,u)7;)] 


(3.51) 


A 

where  the  unit  vector  A:(^,  7)  =  c^x  -b  ajy  +  <^z.  It  does  not  seem  possible,  in  general, 
to  express  the  inverse  time-domain  receiving  characteristic  directly  in  terms  of  the  x  and  y 
components,  and  According  to  (2.21),  for  all  (^,>;)  the  dyadic  receiving 

characteristic  and  its  inverse  satisfy 


-  0  •  Rri^,r,,f)dt'  =  iwHii)  [^(e,7)e(e,»;)  +  ^(e,>7)<^(e,»l)l  (3.52) 


A  A 

where  9{^,r})  and  ^(^,  7)  are  time-independent  unit  vectors  given  by 

^  .  C7 


^Lri)  = 


Cy/^^  +  7* 


and 


k(,’i)  =  - 


”  pi+  ^ 


— — z 


(3.53) 


(3.54) 


VF+F*  "  v/F+F*'' 

Taking  the  inverse  Fourier  transform  of  the  probe-corrected  formula  (2.25)  one  obtains 
the  following  three  alternative  time-domain  probe-corrected  formulas 

1  r+00  r+oo  _ 

7>  0  =  y_^  ■  y_^  y  ^  <  -  «'  +  +  TiyQ)dxodyodt'  (3.55) 
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1  /■  +  <»  /•  +  <»  f  +  OO—^l 

=  7^— ^  /  /  /  Rp  -  i' +  Tiyo)  ■  bp{fQj')dxodyQdt’  (3.56) 

^^.^7r y-oo  J-oo  J-oo 

and 

1  /•  +  <»  /■+00  /•  +  00— _1 

=  77^  /  /  (6 >?.<-<' +  ^a;o  +  J?yo)  ip(ro,<')dt'dio<ij/o  (3.57) 

where  fo  =  lo®  +  yo2/-  As  for  the  acoustic  field,  equation  (3.57)  may  be  the  most  attractive 
for  measurement  purposes  since  the  inner  time  integral  can  be  calculated  at  each  scan  point 
fo  for  the  desired  values  of  Thus,  one  would  not  have  to  store  the  values  of  the 

vector  output  of  the  probe  over  time  to  calculate  the  Radon  transform  of  the  electric  field 
radiated  by  the  test  antenna. 


3.2.1  Formulas  for  Time- Derivative  Probes 


If  the  electromagnetic  probe  is  a  so-called  D-dot  probe  [10],  [11],  [14],  [15]  its  output  from 
an  incoming  plane  wave  is  proportional  to  the  time  derivative  of  the  plane-wave  field  and 
the  frequency-domain  receiving  characteristic  is 

[*<??({.-!)  +  !/<??({,'!)]  (3.58) 

where  the  angularly  dependent  function  Qp{^,  rj)  is  independent  of  the  frequency.  (Note  from 
(2.32)  and  (2.33)  that  a  reciprocal  elementary  electric  dipole  probe  is  a  D-dot  probe.)  With 
this  receiving  characteristic,  the  time-domain  probe-corrected  formula  reduces  to 

=  -^Qp  •  J  ^  J_^  bp{fo,i  +  (xo  +  rtyo)dxodyo  (3.59) 


where  Q.  is  found  from  (2.24)  to  be 


Q~p\^,v) 


X  -  Qp(6?)y] 


HiyV)  =  c^x  +  criy  +  c(z. 


(3.60) 


Combining  the  expression  (3.59)  for  the  time  derivative  of  the  Radon  transform  of  the  near¬ 
field  with  the  expression  (3.47)  for  the  far-field  pattern  shows  that  the  far-field  pattern  is 
given  by 


cos  0  =— 1  r+oo  r+oo  _ 

j^i0,<i>,t)  =  — — Q-  {c~^  cos<i>sin0,c~^  sm<i>sm0)  ■  I  I  bp{fo,t  +  r  ■  fQ/c)dxodyo 

Zirc  ^  J^oo  */-oo 


+  0O  /  +  00  . 


(3.61) 

where  0  <  9  <  ir/2.  Thus,  for  D-dot  probes  the  probe- corrected  formula  simplifies  so  that 
only  the  calculation  of  a  double  integral  is  needed  to  find  the  Radon  transform  T  or  the 
far-field  pattern  T. 
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3.2.2  Formulas  for  Reciprocal  Probes 

Finally  we  consider  the  reciprocal  probe  and  derive  time- domain  reciprocity  relations  between 
the  transmitting  and  receiving  characteristics  of  the  probe,  using  the  Fourier  transform 
along  with  frequency-domain  reciprocity  relations.  We  define  the  time-domain  transmitting 
characteristic  of  the  probe,  when  its  reference  line  is  parallel  to  the  x  axis,  to  be 

the  Ra.don  transform  (3.40)  of  the  field  radiated  by  the  probe  into  the  half  space  z  <  0  when 
it  is  located  at  fo  =  0  and  is  fed  by  the  delta-function  input  signal  ap{i)  =  2irS{t). 

The  frequency-domain  reciprocity  relation  (2.27)  shows  that  for  all  real  u; 


(3.62) 


which  by  use  of  the  inverse  Fourier  transform  gives  the  time-domain  reciprocity  relation 


(3.63) 


where  Ti  is  the  Hilbert  transform  given  in  (3.37).  A  similar  result  holds  when  the  reference 
line  of  the  probe  is  parallel  to  the  y  axis.  Inserting  the  expression  (3.63)  for  and  the  corre¬ 
sponding  expression  for  fVj,  into  (3.51)  gives  a  formula  for  the  inverse  receiving  characteristic 
in  terms  of  the  transmitting  characteristics  in  the  two  orientations. 


3.3  Time-Domain  Sampling  Theorem  and  Numerical 
Far-Field  Calculations 

In  this  section  we  will  consider  only  the  special  cases  where  the  output  of  the  probe,  due 
to  an  incoming  plane  wave,  is  proportional  to  the  time  derivative  of  that  plane-wave  field. 
These  cases  are  of  practical  interest  because  such  probes  au:tually  are  used  for  time-domain 
measurements  [10],  [11],  [14],  [15].  Two  different  computation  schemes,  each  of  which  nu¬ 
merically  calculates  the  time-domain  far-field  pattern  from  sampled  near-field  data,  will  be 
presented.  The  first  is  the  frequency-domain  computation  scheme,  which  is  based  on  the 
frequency-domain  probe-corrected  formulas  derived  in  Chapter  2.  The  second  is  the  time- 
domain  computation  scheme,  which  is  based  on  the  time-domain  probe-corrected  formulas 
derived  in  Chapter  3. 


3.3.1  Frequency-Domain  Computation  Scheme 

This  computation  scheme  is  similar  to  the  non- probe-corrected  scheme  presented  in  [4, 
sec.4.1]  and  consists  of  the  following  three  steps:  (1)  use  the  Fourier  transform  to  calcu¬ 
late  the  frequency-domain  output  of  the  probe  from  the  measured  time-domain  output  of 
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the  probe,  (2)  calculate  the  frequency-domain  far  field  from  the  frequency-domain  output 
of  the  probe,  and  (3)  use  the  inverse  Fourier  transform  to  calculate  the  time-domain  far 
field  from  the  frequency-domain  far  field.  This  scheme  makes  use  of  well-known  frequency- 
domain  probe-corrected  formulas  (2.8)  or  (2.25),  sampling  theorems,  and  the  fast  Fourier 
transform  (FFT).  The  formulas  for  the  discrete  version  of  this  scheme  to  calculate  the  time- 
domain  acoustic  far-field  pattern  T  are  easily  found  from  the  non-probe-corrected  formulas 
[4,  sec.4.1,  eqs.(4.5)-(4.7)].  We  summarize  the  probe-corrected  formulas  for  the  acoustic  field 
below. 

Assume  that  the  output  of  the  probe  is  bandlimited  so  that  can  be  set  equal  to  zero 
for  |cj|  >  uimax  and  assume  that  bp(fQ,t)  begins  at  some  time  to  (which  may  depend  on  the 
position  fo  in  the  scan  plane)  and  ends  approximately  at  a  time  that  can  be  expressed  as 
to  +  —  l)At  where  At  =  7r/(j„,ar.  Then  the  frequency-domain  output  of  the  probe  can 

be  calculated  from 

i>pc^(ro)  =  ^  £  6p(ro,to  +  mAt)e‘"(‘°+"‘^‘>At,  |a;|  <  (3.64) 

fTisO 


The  frequency-domain  far-field  pattern  can  be  found  from  the  probe- corrected  formulas  (2.8) 
and  (2.5).  With  the  receiving  characteristic  given  by  (3.16)  the  discrete  form  of  this  formula 
becomes 


<t>)  =  4)  = 


cos  6 

2ircQp{c~^  cos  4  sin  0,  c"’  sin  ^  sin  0) 

N,  Ny 

•  £  6pw(romn)e*“^''®’"“^"AxoAyo 

m——Ni  nss— 


(3.65) 


where  foron  =  mAxoX  -f-  nAyoy  is  a  sampling  point  on  the  scan  plane.  The  sample  spacing 
Axo  =  Ayo  =  =  ire/ Umax  IS  required  by  the  frequency-domain  sam.pling  theorem 

[1,  fig.  10],  while  Nx  and  Ny  are  integers  that  determine  the  size  of  the  scan  plane,  and 
f  =  X  cos  ^ sin  ^  -I-  y  sin  ^  sin  0  -1-  z  cos  0  is  the  direction  to  the  far-field  observation  point.  Now 
that  the  frequency-domain  far-field  pattern  is  calculated,  the  time-domain  far-field  pattern 
is  found  from  the  inverse  Fourier  transform 


m  4,t)=  E  4)e-'”'^-^Au, 

m=— Af«/2 


Au;  = 


2u„ 


(3.66) 


by  means  of  the  sampling  theorem  if  the  far  field  ^{9,<f>,t)  has  about  the  same  duration  as 
the  output  of  the  probe.  If  the  time  duration  of  the  far  field  is  longer  than  that  of  the  output 
of  the  probe,  one  must  decrease  the  frequency  sample  spacing  Au  by  increasing  The 
numerical  example  given  at  the  end  of  this  section  shows  that  the  duration  of  the  calculated 
far-field  pattern,  in  general,  depends  not  only  on  the  near  field  but  also  on  the  size  of  the 
scan  plane.  This  is  because  the  artificial  edges  of  the  scan  plane  produce  a  diffracted  field 
that  will  make  the  time  duration  of  the  far-field  pattern  calculated  from  (3.66)  longer  than 
that  of  the  exact  far-field  pattern. 
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3.3.2  Time-Domain  Computation  Scheme 

This  computation  scheme  is  similar  to  the  non-probe-corrected  scheme  presented  in  [4, 
sec.4.2]  and  uses  the  direct  formulas  (3,18)  or  (3.61),  which  give  the  time-domain  far-field 
pattern  directly  in  terms  of  the  time-domain  output  of  the  probe.  The  discrete  version  of 
the  formula  (3.18)  to  calculate  the  acoustic  far-field  pattern  J-  is  similar  to  the  non-probe- 
corrected  formula  [4,  sec.4.2,  eq.(4.11)],  and  can  be  found  by  Fourier  transforming  (3.65)  to 
get 

)  27rc(5p(c”*  cos  ^sin  c~*  sin  ^sin  0) 

•  51  6p(fom„,t  +  r  •  fo„in/c)AioAyo-  (3.67) 

m——Nx  n=—Ny 


3.3.3  Comparisons  of  the  Two  Computation  Schemes 

The  comparison  given  in  [4,  sec.4.4]  of  the  non-probe-corrected  computation  schemes  also 
applies  to  the  probe-corrected  ones.  Spf  fifically,  when  the  far  fields  of  electrically  large  radi¬ 
ators  are  calculated  for  all  times  and  all  angles  of  observation,  the  FFT  makes  the  frequency- 
domain  computation  scheme  much  faster  than  the  time-domain  computation  scheme.  When 
only  part  of  the  far  field  is  calculated,  the  difference  in  computer  time  for  the  two  compu¬ 
tation  schemes  becomes  smaller  and  the  time-domain  computation  scheme  becomes  more 
advantageous  because  of  its  simplicity.  The  time-domain  computation  scheme  is  much  eas¬ 
ier  to  use  than  the  frequency-domain  computation  scheme  and  it  is  consequently  the  more 
attractive  scheme  when  one  is  not  concerned  with  the  amount  of  computer  time  it  takes 
to  perform  the  far-field  calculations.  Furthermore,  the  calculated  far-field  pattern  has  an 
extended  time  duration  because  of  the  finite  scan  plane  and  to  avoid  aliasing,  this  extended 
duration  has  to  be  taken  into  account  when  the  frequency-domain  computation  scheme  is 
used.  The  time-domain  computation  scheme  does  not  make  use  of  the  Fourier  transform 
and  thus  does  not  encounter  this  aliasing  problem. 

As  pointed  out  in  Section  3.1,  planar  time-domain  near- field  antenna  measurements, 
unlike  single-frequency  near-field  measurements,  have  the  capability  of  eliminating  finite 
scan  errors.  This  makes  it  possible  to  determine  the  far-field  pattern  of  broadbeam  antennas 
from  planar  time-domain  measurements,  as  will  be  demonstrated  in  the  following  example. 


3.3.4  Numerical  Example 

The  probe-corrected  formula  (3.18),  and  its  discrete  version  (3.67),  will  now  be  verified 
numerically  by  using  the  time-domain  computation  scheme  to  calculate  the  far-field  pattern 
of  an  acoustic  point  source  with  Gaussian  time  dependence.  The  near-field  data  used  in  this 
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computation  are  obtained  by  measuring  the  near  field  of  the  point  source  with  a  nonidea] 
probe. 

The  acoustic  point  source  with  Gaussian  time  dependence  is  located  at  fj  =  ~dz  and  its 
radiated  field  is  given  by 


fit  -  +  d^l/c) 

47r|r  +  dz\  ’ 


m  = 


(3,68) 


where  t  and  d  are  positive  constants.  The  near-field  data  are  obtained  by  using  a  probe 
whose  receiving  characteristic  is  in  the  form  (3.31),  that  is,  the  output  of  the  probe  due  to  the 
incoming  plane  wave  is  proportional  to  the  time  derivative  of  the  incoming  plane  wave 
The  angularly  dependent  function  Qpi(,T))  is  chosen  to  be  cos  6,  where  0  is  the  angle 
between  the  propagation  direction  of  the  plane  wave  and  the  z  axis.  Since  the  propagation 
direction  of  the  plane  wave  is  given  hy  k  =  ^cx  +  rjcy  +  ^cz,  the  angularly  dependent  function 
is 


Qpi(.v)  =  <  = 


( 


lv/1  +T)^  <C  2 

—  1|, 


(3.69) 


To  calculate  the  double  summation  (3.67)  for  the  far-field  pattern,  we  have  to  know  the 
output  of  the  probe  at  each  scan  point  on  the  scan  plane.  Since  the  receiving  characteristic 
of  the  probe  is  defined  in  terms  of  plane  waves  we  have  to  expand  the  acoustic  point-source 
field  in  terms  of  plane  waves  in  order  to  determine  the  output  of  the  probe. 

We  start  by  writing  the  acoustic  point-source  field  (3.68)  in  terms  of  the  Fourier  integral 


+0O 

“OO  47r|r  -I-  dz\ 


(3.70) 


where  =  r(4>/7r)~’e~‘*’^’‘*/^®  is  the  spectrum  of  the  Gaussian  time  function  /(t)  in  (3.68). 
The  reason  we  write  the  acoustic  field  in  this  form  is  that  we  can  now  make  use  of  the 
identity  [4,  eq.(2.29)]  to  show  that 


47r|f  -|-  dz( 


Stt'  J-oo  J-oo  7 


(3.71) 


which  gives  the  frequency- domain  point-source  field  in  terms  of  plane  waves.  The  plane-wave 
spectrum  for  the  point-source  field  (3.70)  is  then  seen  from  (3.71)  and  (2.1)  to  be 


To^ik 


xt  ^y) 


(3.72) 


Then  (2.7)  shows  that  the  frequency-domain  output  of  the  probe  is  given  by 

bpu>iro)  =  ^  r°°  ^  0  (3.73) 

47r  J—oo  J—oo  7 
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where  Rp^  is  the  frequency- domain  receiving  characteristic  of  the  special  time-derivative 
probe  found  from  (3.16)  and  (3.69)  to  be 


Rpti!  (  )  — 


2ir 


27r 


ic-y 

2^' 


(3.74) 


Inserting  this  expression  for  the  receiving  characteristic  of  the  probe  into  (3.73)  and  making 
use  of  (3.71),  one  finds  that  the  frequency-domain  output  of  the  probe  is  given  by 


rf  8  f,ik\fo-{-di\ 
“  “4x  dd\fa  +  dz\' 


(3.75) 


Note  that  this  expression  could  also  have  been  given  in  terms  of  a  partial  derivative  with 
respect  to  z.  However,  since  ^  =  0  on  the  scan  plane  it  is  more  convenient  to  write  it  in 
terms  of  the  distance  d  between  the  point  source  and  the  scan  plane.  Taking  the  inverse 
Fourier  transform  of  (3.75)  shows  that  the  time-domain  output  of  the  probe  is  simply 

(3-76) 

where  is  the  Gaussian  acoustic  point-source  field  (3.68).  We  have  now  shown  that  the 
output  of  the  probe  equals  — c  times  the  partial  derivative  with  respect  to  d  (or  z)  of  the 
point-source  field.  This  result  could  also  have  been  obtained  by  using  the  formulas  of  [7] 
where  receiving  antennas  are  described  as  linear  differential  operators.  When  the  distance 
between  the  probe  and  the  point  source  is  large,  (3.76)  shows  that  the  output  of  the  probe  is 
6p(ro.  t)  ^  cos  t)  where  9  is  the  angle  between  the  z  axis  and  the  vector  going  from 

the  point  source  to  the  point  fo  on  the  scan  plane.  This  result  is  explained  by  noting  that 
the  point-source  field  at  a  position  fo  in  the  far  field  consists  of  a  plane  wave  propagating 
in  the  direction  given  by  9,  and  the  output  of  the  special  probe  equals  cos  9  times  the  time 
derivative  of  an  incident  plane-wave  field. 

Having  calculated  the  output  of  the  probe  we  use  (3.18)  to  find  the  following  expression 
for  the  far-field  pattern  of  the  test  antenna 

J='{9,  (t>,  t)  =  /  /  — $®(ro,<  +  r  •  rofc)dxodyo,  0<9  <  ir/2.  (3.77) 

The  discrete  version  of  (3.77)  useful  for  computation  is  obtained  from  (3.67)  as 
1  /V,  B 

J^{9,4>,t)  =  52  52  +  r  ■  fomn/c)AxoAyo,  0  <9  <  7r/2.  (3.78) 

m=-N^  n=-Ny 

If  probe  correction  had  been  neglected,  the  far-field  pattern  would  have  been  given  by 
the  right  side  of  (3.77)  or  (3.78)  divided  by  cos^.  Thus,  probe  correction  can  be  neglected 
only  if  cos  9  1. 

Let  us  now  show  the  computed  far-field  pattern  of  the  Gaussian  acoustic  point  source 
obtained  from  the  near-field  formula  (3.78).  In  the  previous  report  [4,  sec.4.3]  it  was  shown 
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Figure  3.5:  Gaussian  point  source  measured  on  a  finite  scan  plane  with  a  time-derivative 
probe. 

that  the  maximum  effective  frequency  in  the  acoustic  field  (3.68)  is  u^max  —  12/t  and  thus 
the  shortest  wave  length  is  =  2irc/u;„ax  —  ctI2.  For  the  numerical  calculations  we 
choose  d  =  ~  cr.  To  satisfy  the  requirements  of  the  sampling  theorem  [4,  sec. 4. 2],  the 

near  field  is  sampled  with  spatial  sample  spacing  Axo  =  Aj/o  =  Ami„/2  ~  ct/4.  The  far- field 
pattern  will  be  calculated  for  —  (45°,  0)  and  the  scan  plane  is  chosen  to  be  a  square 

of  sidelength  20d,  which  is  large  enough  to  avoid  the  problems  of  the  interference  between 
the  correct  field  and  the  erroneous  field  for  this  angle  of  observation.  The  erroneous  field  is 
due  to  the  edges  of  the  scan  plane  and  has  been  discussed  in  (4,  sec.4.4].  With  this  size  of 
the  scan  plane  and  this  observation  direction,  the  erroneous  signal  ends  approximately  20. 5r 
after  the  correct  signal.  Thus,  the  duration  of  the  far-field  pattern  calculated  from  (3.77) 
is  approximately  23r  whereas  the  duration  of  the  exact  far-field  pattern  is  only  2.5t.  See 
Figure  3.5  for  a  picture  of  the  scanning  geometry. 

Figure  3.6  shows  the  far-field  pattern  calculated  for  (9,<i>)  =  (45°,  0)  in  three  different 
ways:  (a)  is  the  exact  Gaussian  far-field  pattern;  (b)  is  the  far-field  pattern  calculated  without 
probe  correction  by  inserting  the  output  of  the  probe  (3.76)  as  if  it  were  the  measured  near 
field  into  [4,  eq.(4.11)];  and  (c)  is  the  far-field  pattern  calculated  by  the  probe-corrected 
formula  (3.78). 

The  curve  (a),  which  is  exact,  has  the  Gaussian  wave  form  and  is  only  significantly  nonzero 
on  the  interval  — 0.5r  <  t  <  2. Or.  The  curve  (b),  which  is  obtained  without  probe  correction. 
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Figure  3.6:  Values  of  the  far-field  pattern  .^(45°,  0,<).  (a)  exact;  (b)  without  probe  correction; 
(c)  with  probe  correction. 

is  clearly  in  error  even  on  the  interval  where  the  exact  far-field  pattern  is  nonzero.  In  fact,  it 
equals  cos 45°  «  0.707  times  the  probe-corrected  value  given  by  the  curve  (c).  The  curve  (c), 
which  is  obtained  from  the  probe-corrected  formula  (3.78),  cannot  be  distinguished  from  the 
exact  curve  on  the  interval  — r  <  t  <  2.4t.  Because  of  the  finite  size  of  the  scan  plane  the 
curves  (b)  and  (c)  are  erroneous  (negative)  on  parts  of  the  interval  from  t  =  2.4t  to  t  =  23t 
(the  end  of  the  erroneous  signal  {t  =  23t)  is  found  by  determining  the  time  after  which  the 
field  is  identically  zero  on  the  finite  scan  plane  [4,  sec.4.4)).  However,  the  erroneous  parts 
for  the  curve  (c)  do  not  overlap  with  the  interval  where  the  exact  far-field  pattern  is  nonzero 
so  the  curve  (c)  is  an  excellent  approximation  to  the  exact  far-field  pattern.  The  curve  (b) 
obtained  without  probe-correction  is  seen  to  be  a  very  poor  approximation  to  the  exact  far- 
field  pattern.  Probe  correction  is  therefore  required  to  get  an  accurate  far-field  pattern  when 
this  special  time-derivative  probe  with  cos  9  angular  dependence  is  used.  Furthermore,  it  is 
demonstrated  that  time-domain  planar  near-field  measurements  can  be  used  to  accurately 
determine  the  far-field  pattern  of  broadbeam  antennas  that  cannot  be  accurately  measured 
in  the  frequency  domain  because  of  the  errors  introduced  by  the  finite  scan  plane. 
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OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  ReNabiiity  Sciences. 


